Cold Bose gases with large scattering lengths 
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We calculate the energy and condensate fraction for a dense system of bosons interacting through 
an attractive short range interaction with positive s-wave scattering length a. At high densities, 
n 3> a~ 3 , the energy per particle, chemical potential, and square of the sound speed are indepen- 
dent of the scattering length and proportional to n 2//3 , as in Fermi systems. The condensate is 
quenched at densities na 3 ~ 1. 
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The densities of atomic vapors under experimental con- 
ditions are generally low, in the sense that the distance 
between particles is typically much larger than the mag- 
nitude of the atomic scattering length a, which character- 
izes the strength of two-body interactions. Recently, it 
has become possible to tune atomic scattering lengths to 
essentially any value, positive or negative, by exploiting 
Feshbach resonances Q . This opens up the possibility of 
creating systems which are dense, in the sense that the 
"gaseousness parameter", n|a| 3 , where n is the particle 
density, is large compared with unity j| . 

In dilute alkali gases, the bare two-body interaction 
at distances large compared with the atomic size is at- 
tractive due to the van der Waals interaction, but the 
effective two-body interaction at low energies, which is 
given by 4irh 2 a/m, may be large and positive due to the 
existence of a bound state just below threshold. 

One of the novel features of such gases is that the char- 
acteristic length R associated with the range of the in- 
teratomic potential is of order 10 -6 cm. The quantity 
ro = (3/47T71) 1 / 3 , which is a measure of the atomic sepa- 
ration, is typically of order 10 -4 cm, and therefore much 
larger than R. In this Letter we consider the properties 
of a Bose gas for arbitrary n\a\ 3 , assuming that R <C r Q . 
The corresponding problem for a Fermi gas has been con- 
sidered previously H. 

We begin by reviewing briefly the properties of the 
dilute Bose gas with repulsive interactions, a > 0. At 
low densities (na 3 < 1) the interaction energy is propor- 
tional to the scattering length and the density, as was 
first found by Lenz Q. With higher order terms calcu- 
lated by Huang, Lee, and Yang H and Wu Q included, 
the energy per particle is 
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To this order, the only property of the two-body inter- 
action that enters is the scattering length. These results 



have been used to calculate the leading corrections be- 
yond mean-field theory to the properties of trapped Bose 
gases . Terms of order na 3 and higher depend on other 
properties of the interaction, as discussed in detail re- 
cently in Ref. [|. 

The expansions leading to the above results are asymp- 
totic, and consequently they give little guidance to the 
properties of a dense gas, na 3 ^l. Even though the range 
of the interaction is small compared with the particle 
separation, correlations between particles are very im- 
portant, as one can see from the fact that all terms in 
the low-density expansion Eq. (uh are formally of the 
same order of magnitude when nar ~ 1. 

The Jastrow wave function [Rill 



r(ri, ...,rjv) = Y[f(i 
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provides a good approximation for cold dense Bose sys- 
tems such as the hard-sphere gas and liquid 4 He. In 
these systems the pair correlation function f(r) may be 
determined variationally by minimizing the expectation 
value of the energy, E/N = (^\H\^) / which 
may be calculated by Monte Carlo methods pl[ . The 
hard-sphere gas problem has also been solved, essentially 
exactly, by quantum Monte Carlo methods for scat- 
tering lengths almost as large as ro. 

A much simpler method, the lowest order constrained 
variational (LOCV) method Jl3| ], provides a fair approx- 
imation to the results of quantum Monte Carlo cal- 
culations. The basic idea of this method is that for 
r < ro the Jastrow function /(r) approximately obeys 
the Schrodinger equation for a pair of particles, 
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m dr 



rf(r) = Xrf(r) , 



(3) 



where A is the pair energy, while for r 3> ^o, / tends 
to unity. To take into account many-body effects, which 
become important when r is comparable to ro, we impose 
the condition that f(r > d) = 1 and f'(r = d) = 
at some healing distance d. In the LOCV method the 
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expectation value of the energy is approximated by the 
two-body cluster contribution 



E/N = 27mA / f 2 (r)r 2 dr 
Jo 



and d is chosen such that 
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f 2 (r)r 2 dr = 1 
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This ensures that on average there are only correlated 
pairs, and E/N = A/2. The energy of a pair is A, and 
thus the calculation of the energy by evaluating the ex- 
pectation value of the Hamiltonian is consistent with the 
eigenvalue of the Schrodinger equation (||) ■ 

When the range of the interaction is small compared 
to both \a\ and d, as in dilute alkali gases, the interaction 
can be replaced by the boundary condition (rf)'/rf = 
— 1/a at r = 0. This boundary condition is also satisfied 
by the two-body s-wave scattering wave function, and 
the scattering length a is then the position of the node 
of f(r) at A = 0. The solution to Eq. (||) for positive 
A = k 2 /to, is rf(r) oc sin[fc(r — b)], with a node at r = b, 
and kb = arctanfca. In the dilute limit (ka — » 0), b is 
equal to a, and in the dense limit (ka — > oo), b — n/2k. 

The other boundary condition, f'(r = d) = 0, gives 



t tan k — 1 
1 + k tan k 
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where n = kd, and the value of d is determined from 
Eq. (g). When / 2 (r) = 1, Eq. (g) gives d = r . Since 
f(r ru d) ~ 1 the value of d is generally close to ro. 

In the LOCV approximation /(r) for a short range 
interaction is therefore given by 



f(r < d) 



d sin[fc(i — b)] 
r sin[fc(c2 — b)] 
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It is negative for r < b, and changes sign at r = b. The 
Jastrow wave function (Eq. ||) does not describe the node- 
less ground state in this case. When a <C r the gas has 
only pairwise collisions, and the Jastrow wave function 
^fj describes its state accurately. We assume that, as 
the gas is compressed or a is increased adiabatically, the 
nodal structure of the wave function is maintained and 
that \& j continues to be a good approximation to the true 
wave function. However, as a or n increases, the rate of 
three-body collisions, which deplete the monoatomic gas 
density by forming molecules, will increase and limit the 
lifetime of the monoatomic gas. 

The variational principle cannot be used to determine 
the Jastrow function f(r) in the case of attractive short 
range interactions, because the gas is not in its ground 
state. In this case, it may easily be verified that the 
energy of the g function of the healing distance d, 

has an uninteresting minimum at d —> with value zero. 
Such a minimum exists for all ^-function interactions. 



However, we anticipate that for the present problem the 
LOCV method will provide a good approximation if we 
take Eq. (^) for the Jastrow factor. 

In Fig.l we show the energy per particle for hard- 
sphereand attractive short range potentials in the LOCV 
approximation, along with the low-density expansion (|l|) 
(without the logarithmic term) and the Lenz expression 
as a function of a/ro, and the exact quantum Monte Carlo 
results for hard-spheres. The LOCV energy is close to the 
low density expansion at small values of a/ro, though it 
does not have the correct analytical form. It contains 
the leading Lenz term in Eq. (nh, as may be verified by 
expanding the tan k in Eq. ^ in powers of k |14[ , but the 
(na 3 ) 1 / 2 and higher terms are approximated. The LOCV 
energies for the hard-sphere gas are in excellent agree- 
ment with the exact ones for a/ro < 0.6, but for larger 
a/ro the results for the hard-sphere and the attractive 
short range potentials deviate significantly. 




FIG. 1. Energy per particle, multiplied by 2mrg/ft 2 vs. 
a/ro. The LOCV results for attractive short range and 
hard-sphere interactions are shown by full and dashed lines re- 
spectively, while the dash-dot and dotted lines show the Lenz 
result and the low density expansion. Monte-Carlo numerical 
results for a hard-sphere gas are shown as stars. 



In the small a/ro limit, using the result b = a, and 
d ~ ro we find that the LOCV correlation function is 
given by 



f(r < r Q ) 



ro sin[fc(r — a)] 
r sin[fc(ro — a)] 
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with k 2 = 3a/rQ. At small r this correlation is ~ 1 — a/r, 
as in the two-body scattering state at zero energy. 

By contrast, when a/ro is large we have b — n/2k and 
Eq. (H) reduces to 



k tan k = — 1 , 



(9) 
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with solutions kd = n = 2.798386.. 
asymptotically k — mi for integer n. 
lation function is given by 

d cos kr 

fir <d) = — . 

r cos kd 



6.1212..,..., and 

The pair corre- 



(10) 



Substituting this expression in Eq. (g) gives d = 
(2/3) 1 /'V for kd = 2.798386..., appropriate for a Jas- 
trow factor / with one node. The energy is given by 
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This energy is proportional to n 2 ^ as one would expect 
on dimensional arguments. A similar result, (E/N oc 
n 2 / 3 ), was found for fermions B in the a — > — oo limit. 
(We obtain E/N = lb.26h 2 n 2 >^/m on approximating d 
by r .) 

For the case of attractive short range interactions it 
is difficult to go beyond the LOCV calculation. When 
a/ro < 0.2 we can perform Monte Carlo calculations 
of the expectation value of the energy, and the results 
are not significantly different from those from the LOCV 
method. At larger values of a/ro the state described by 
the Vtj has a clustering instability. All the atoms end up 
in a few dense clusters as the Metropolis walk progresses, 
as is expected to happen with time in experiments. 

The rate of molecule formation is proportional to the 
probability of finding three atoms together. In the lowest 
order (3-body) cluster approximation the probability of 
finding three atoms at points rj, r 2 , and r 3 is propor- 
tional to / 2 (ri2)/ 2 (r 2 3)/ 2 (r 31 ), where r l3 = |r 4 - rd. It 
diverges as all — > for the / given by Eq. (|10) or 
Eq. (^). However, within the interaction radius R the 
true / does not have the 1/r behavior. We can estimate 
the probability P 3 (i?) of finding two atoms within the 
interaction range R of a third by approximating the / 
at r < R by f(R). In the dense limit Eq. (10) gives 
f(r<R)~ —ro/R, and 
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independent of the gas density. In contrast, at small 
values of a/r Eq. (m gives f(r < R) ~ — a/ro, an d 
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At present it is unclear how the behavior can be matched 
on to the a 4 dependence predicted for n\a\ 3 <C 1 (see, e.g., 
Ref. @). 

We next estimate the depletion of the condensate. 
In the LOCV approximation, the condensate fraction is 
given by @ 
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[l-f(r)] 2 d 3 r. 



(14) 



Substituting Eq. (^) in the above equation gives 
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— (sin k — k cos k) — 1 



(15) 



The condensate fraction calculated using the LOCV ap- 
proximation is shown in Fig. 2 along with the predictions 
for the dilute limit given by pi 
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Results for the hard-sphere potential obtained with the 
LOCV method, and exact quantum Monte Carlo calcu- 
lations jl^] are also shown in Fig. 2. The LOCV results 
for no/n are rather crude, since long-range correlations 
neglected in the LOCV approximation have a significant 
effect on no. All the results indicate that superfluidity 
will be quenched when a^ro- In the case of hard-spheres 
it is known that the liquid solidifies at ro ~ a |L5j . 
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FIG. 2. Condensate fraction vs. a/ro. LOCV calculations 
for attractive short range (full line) and hard-sphere (dashed 
line) interactions, Eq. [L6j (dotted line) and exact QMC results 
for hard-spheres (stars). 



The chemical potential at zero temperature is given by 
fi = dE/dN . In the dilute limit the leading contributions 
to (i are \x = Airh 2 na/m [1 + (32/3)(na 3 /7r) 1 / 2 ], while, in 
the dense limit, the LOCV approximation gives 



(j, = 22.22- 



h 2 n 2 ' 3 



(17) 



The calculated fi(n) is shown in Fig. 3. 

The square of the sound speed c s is given by c 2 = 
(n/m)d[i/dn, and this is also shown in Fig. 3. Its di- 
lute and dense limits are given by c 2 = Anf^na/m 2 [1 + 
^na 3 /?!-) 1 / 2 ] and 
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14.81 ■ 



Ti 2 n 2 ' 3 



(18) 



respectively. At high densities the collective modes in 
the hydrodynamic limit are similar to those for fermions 
fl8f since the energy per particle scales with density 
as n 2 / 3 , as it does for a free Fermi gas. The energy 
per particle in a one-component Fermi gas, E/N = 
(3/10)h 2 (6ir 2 n) 2 / 3 /m, is smaller by a factor 2.92 than 
the LOCV result in the dense Bose limit. The chemical 
potential and sound speed squared are therefore corre- 
spondingly larger in a dense Bose gas. 
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FIG. 3. The zero temperature chemical potential mul- 
tiplied by 2mro/?i 2 for: attractive short range (solid line) 
and hard-sphere (dashed line) from LOCV calculations, low 
density expansion (dotted line) and the square of hydro- 
dynamic sound speed multiplied by 2m 2 ro/7i 2 for: LOCV 
for attractive short range (double dash-dot), LOCV for 
hard-sphere(dash-double dot) and low density expansion 
(dash-dot). 



The calculations of the energy density and the estimate 
of the three-body cluster probability P3 (R) indicate that 
in limit where the scattering length is large compared to 
the interparticle spacing, the scattering length becomes 
irrelevant. These predictions could be investigated exper- 
imentally by measuring the size of a cloud of condensate 
and the particle lifetime as the scattering length is varied 
by exploiting a Feshbach resonance. 
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